Longitudinal double-spin asymmetry and cross section for inclusive neutral pion 
production at midrapidity in polarized proton collisions at ^/s = 200 GeV 
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We report a measurement of the longitudinal double-spin asymmetry All and the differential cross 
section for inclusive tt production at midrapidity in polarized proton collisions at ^fs = 200 GeV. 
The cross section was measured over a transverse momentum range of 1 < pr < 17 GeV/c and 
found to be in good agreement with a next-to-leading order perturbative QCD calculation. The 
longitudinal double-spin asymmetry was measured in the range of 3.7 < pr < 11 GeV/c and excludes 
a maximal positive gluon polarization in the proton. The mean transverse momentum fraction of 
7T°'s in their parent jets was found to be around 0.7 for electromagnetically triggered events. 

PACS numbers: 13.87.Ce, 13.88. +e, 14.70.Dj, 12.38.Qk 
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The spin structure of the nucleon is one of the fundamen- 
tal and unresolved questions in Quantum Chromodynam- 
ics (QCD). Deep-inelastic scattering (DIS) experiments 
studying polarized leptons scattered off polarized nuclei 
have found the quark and anti-quark spin contributions 
to the overall spin of the nucleon to be small, at the level 
of 25% [TJ |2] , leading to increased interest in the spin 
contribution from gluons. DIS experiments have placed 
coarse constraints on the polarized gluon distribution func- 
tion Ag(x), based on the scale dependence of polarized 
structure functions [3] 0] and on recent semi-inclusive 
data [3 |7]. Measurements using collisions of longi- 
tudinally polarized protons are attractive because they 
provide sensitivity to the polarized gluon spin distribution 
at leading order through quark-gluon and gluon-gluon 
scattering contributions to the cross section. 

The sensitivity of inclusive hadron and jet produc- 
tion to the underlying gluon polarization in high-energy 
polarized proton collisions has been discussed in detail 
in Refs. [5J [9]. The theoretical framework in the con- 
text of next-to-leading order perturbative QCD (NLO 
pQCD) calculations is very well developed to constrain 
Ag(x). The first global analysis of semi-inclusive and 
inclusive DIS data, as well as results obtained by the 
PHENIX PJ]J and STAR [11J experiments, placed a 
strong constraint on Ag{x) in the gluon momentum- 
fraction range of 0.05 < x < 0.2, and suggested that 
the gluon spin contribution is not large in that range [T2"] . 
This conclusion was driven primarily by data on inclusive 
hadron and jet production in polarized proton collisions 
at y/s = 200 GeV at RHIC. 

In this paper, we report on the measurement of the 
cross section and the longitudinal double-spin asymmetry 
All for inclusive 7r° production at midrapidity in polar- 
ized proton collisions at s/s = 200 GeV by the STAR 
experiment [12] at RHIC. The cross section is compared 
to a NLO pQCD calculation and the observed agree- 
ment provides an important basis to apply pQCD for the 
interpretation of All- The asymmetry is defined as 



where a ++ and c H are the inclusive 7r° cross sections 
for equal (++) and opposite (H — ) beam helicity configu- 
rations. The measured longitudinal double-spin asymme- 
try probes a gluon momentum fraction of approximately 
0.03 < x < 0.3, and is compared to NLO pQCD cal- 
culations. In addition, we present the mean transverse 
momentum fraction of 7r°'s in electromagnetically trig- 
gered jets. This measurement allows one to relate the 
spin asymmetry measurements performed with inclusive 
7r°'s to those using reconstructed jets. It may also help 
to constrain fragmentation models. 

The data for the analyses presented here were collected 
at STAR in 2005 using stored polarized 100 GeV proton 
beams with an average luminosity of 6 x 10 30 cm" 2 s _1 . 
Longitudinal polarization of proton beams in the STAR 
interaction region (IR) was achieved by spin rotator mag- 
nets upstream and downstream of the IR that changed the 
proton spin orientation from its stable vertical direction to 



longitudinal [Tl]. The helicities were alternated between 
successive proton bunches in one beam and pairs of suc- 
cessive proton bunches in the other beam. This allowed 
us to obtain all four helicity combinations of the colliding 
bunch pairs at the STAR IR in quick succession. Addi- 
tional reduction of systematic uncertainties was achieved 
by periodically changing the helicity patterns of the stored 
beams. The polarization of each beam was measured sev- 
eral times per fill using Coulomb-Nuclear Interference 
(CNI) proton-carbon polarimeters |15j . which were cali- 
brated using a polarized hydrogen gas-jet target [16] . The 
average RHIC beam polarizations in the 2005 run were 
Pi = 52 ± 3% and P 2 = 48 ± 3%. Non-longitudinal beam 
polarization components were continuously monitored 
with local polarimeters at STAR [17] and were found to 
be no larger than 9% in absolute magnitude. 

The principal STAR detector subsystems for the mea- 
surements presented here were the Barrel Electromagnetic 
Calorimeter (BEMC) [TB] and the Beam-Beam Counters 
(BBC) [T7]. In addition, the Time Projection Cham- 
ber (TPC) [13] was used for vertexing, for measuring 
the charged component in the reconstructed jets, and as 
a charged particle veto for the ir° reconstruction. The 
BEMC is a lead-scintillator sampling calorimeter with 
a granularity of At] x Aip — 0.05 x 0.05 rad, where one 
such cell is referred to as a tower. It contains a shower 
maximum detector (BSMD) that consists of two layers 
of wire proportional counters with cathode strip readout, 
one in the azimuthal direction and one in the longitu- 
dinal direction, at a depth of about 5 radiation lengths 
in each calorimeter module, providing a segmentation of 
0.007 x 0.007 rad. For the 2005 running period, half of 
the BEMC was instrumented and operational, provid- 
ing 2ir azimuthal coverage for < r\ < 1. The BBCs 
are composed of segmented scintillator rings, covering 
3.3 < \rj\ < 5.0 on both sides of the IR. The BBCs 
were used to trigger on collisions, to measure the helicity- 
dependent relative luminosities, and to serve as local 
polarimeters. The TPC provided charged particle track- 
ing inside a 0.5 T solenoidal magnetic field over the full 
range of azimuthal angles for |t^| < 1.3. 

Proton-proton collisions in the STAR detector were 
identified by a minimum bias trigger (MB), defined as a 
coincidence of hits in both BBCs. The cross section for 
this trigger was ctbbc = 26.1 ± 0.2 (stat) ± 1.8 (syst) mb, 
corresponding to 87 ±8% of the non-singly diffractive p+p 
cross section at y/s = 200 GeV ^ZU\ . Rare hard scattering 
events were selected by two high-tower triggers, HT1 and 
HT2, that required a transverse energy deposition in a 
single BEMC tower above thresholds of 2.6 and 3.5 GeV, 
respectively, in addition to satisfying the MB condition. 

A data sample with an integrated luminosity of C = 
0.17 nb _1 for MB, 0.16 pb _1 for HT1, and 0.66 pb _1 for 
HT2 triggers was analyzed for the inclusive cross sec- 
tion measurement. Data with an integrated luminosity 
of 0.4(2.0) pb _1 of HT1 (HT2) triggers were used for 
the All determination. The event selection criteria for 
the asymmetry analysis were identical to those used in a 
previously published jet measurement [TT]. About 22% 
of HT1/HT2 triggered events also entered the jet All 
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measurement [TT] , but represented a negligible fraction 
of the much larger inclusive jet data set. Therefore, the 
statistical correlation of the present 7r° and jet All mea- 
surements is negligible. 

Neutral pions were reconstructed in the decay channel 
7r — ► 77 in an invariant mass analysis of pairs of neutral 
BEMC clusters, i.e., those that did not have a TPC track 
pointing to them, with a cut on the two-particle energy 
asymmetry of \E\ — i? 2 |/(£'i+-E2) < 0.7. The tower granu- 
larity was insufficient to resolve cluster pairs in HT1/HT2 
data because of the small opening angle between daughter 
photons of pions that satisfied these triggers. Therefore, 
the BSMD clusters were used to determine the photon 
coordinates in those data. A fiducial volume cut on the 
detector pseudorapidity of 0.1 < r\ < 0.9 was imposed. 
The reconstructed value of the pion pseudorapidity with 
respect to the vertex position was required to fall in the 
range < T] < 1. The 7T° yield was extracted in px bins by 
integrating the background-subtracted invariant mass dis- 
tribution in a ^-dependent window around the tt° peak 
that corresponded to an approximately ± 3 a range. The 
combinatorial background was determined using the event 
mixing method with a jet alignment correction |21l 122] . 

The cross section for ir° production is given by 

F— 1 — 

~ 2ir PT Ap T Ari C £' ( ' 

where Apt and A77 are the bin widths in px and pseudo- 
rapidity, N is the 7T° yield in a bin, and c is an overall 
correction factor that accounts for acceptance, recon- 
struction, and trigger efficiency in that bin, which was 
determined using a Monte Carlo simulation of 7r°'s passed 
through the GEANT [23] model of the STAR detector. 
Figure [T] shows the differential cross section for inclusive 
7r° production. This analysis covered the pion transverse 
momentum range of 1 < px < 17 GeV/c, and data points 
were scaled to the bin centers using local exponential fits 
around each bin. The cross sections up to 4 GeV/c were 
measured using MB triggered events; above 4 (7) GeV/c 
the entries were obtained from HT1 (HT2) triggers. The 
different trigger samples agreed within errors. 

The dominant systematic uncertainty (25% on average) 
of the measured cross section was due to a 5% uncertainty 
in the global energy scale of the BEMC. The other sys- 
tematic uncertainties were related to yield extraction (7%), 
reconstruction efficiency (6%), and relative normalization 
of HT1/HT2 and MB triggers (5%). An additional un- 
certainty due to the limited quality of the electromagnetic 
shower simulation at low photon energies in our GEANT 
model was assigned to the cross section obtained from 
HT1/HT2 data [15(2)% at p T = 4(7) GeV/c]. 

In Fig. [l] the measured cross section is compared 
to a NLO pQCD calculation [8] performed using the 
CTEQ6M set of unpolarized parton distribution func- 
tions [25] and the DSS set of fragmentation functions [24"] . 
In this calculation, the factorization and renormalization 
scales were identified with px (solid curve), and were 
varied by a factor of two to estimate the impact of scale 
uncertainties (dashed curves). The DSS analysis included 
recent measurements of ir° production at midrapidity by 
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FIG. 1: (color online) (a) Cross section for inclusive 7r° pro- 
duction at midrapidity in p + p collisions at y/s = 200 GeV, 
compared to a NLO pQCD calculation [S] based on the DSS 
set of fragmentation functions [24], and to the STAR 71-* mea- 
surement [25] . (b) The ratio of measured cross section and 
the NLO pQCD calculation. The scale uncertainty is indi- 
cated by the dashed curves (p, = 2px, pt/2). The error bars 
are statistical and shaded bands are pr-correlated systematic 
uncertainties. The normalization uncertainty is indicated by 
a shaded band around unity on the right-hand side. 

PHENIX [10J and at forward rapidity by STAR [27]. 
The NLO pQCD calculation shows, within errors, good 
agreement with our data in the fragmentation region 
px > 2 GeV/c. We also compare the cross section for 7r° 
production to the STAR tt measurement [25]. The 7r° 
and (tt + + n~ ) /2 cross sections are expected to be equal, 
and the two STAR measurements agree within statistical 
errors, in spite of using independent detector sub-systems. 

The transverse momentum fraction carried by a high- 
px 7T° in its parent jet, z = px(ir )/pT(jet), was investi- 
gated by associating pions with jets found in the same 
event |28j . The 7r° sample, defined by the invariant mass 
window, contained ss8% of combinatorial background. 
An association was made if the pion was within a cone of 
radius R = J (A?/) 2 + (A<^) 2 = 0.4 around the jet axis. 
The analysis was restricted to 0.4 < r\ < 0.6 in the jet 
pseudorapidity, so that the reconstructed jets were fully 
contained in the BEMC acceptance. The transverse mo- 
mentum of the jet was required to exceed 5 GeV/c. The 
jet was required to have a neutral energy fraction less 
than 0.95, in order to minimize contributions from beam 
background to the reconstructed jet sample. 

Figure |2ja) shows the mean value of z as a function 
of pion px, combined for HT1 and HT2 triggers. The 
data points are plotted at the bin centers in pion px- The 
results were not corrected for detector effects, such as ac- 
ceptance, efficiency, or resolution of the jet reconstruction. 
The systematic error band shown includes contributions 
from the uncertainty of the jet energy scale, the influence 
of the cut on minimum jet px, the contribution of events 
with z > 1, and a variation of other analysis cuts. 

The (z) of 7r°'s in electromagnetically triggered jets was 
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FIG. 2: (color online) (a) Mean transverse momentum fraction 
of 7r°'s in their associated jets, as a function of pion pr, for 
electromagnetically triggered events. Systematic errors are 
shown by the shaded band around the data points. The curves 
are results from simulations with the pythia event generator. 
The solid curve includes detector effects simulated by GEANT, 
while the dashed curve uses jet finding at the pythia particle 
level, (b) The distribution of z for one pr bin, compared to 
pythia with a full detector response simulation. 



found to be around 0.7 and to rise slightly with pion pr, 
consistent with measurements of leading charged hadrons 
in jets in fixed-target experiments The results also 
compare well to recent theoretical calculations for charged 
pions [30] , considering the increase of the measured pion 
momentum fraction due to energy not reconstructed in 
the jet. The expectations from a PYTHiA-based (version 
6.205 [2] with 'CDF Tune A' settings [32]) Monte Carlo 
simulation are also shown. The (z) measured in jets found 
on the pythia particle level, i.e., without any detector 
effects, is lower than in the data due to resolution effects 
and losses in the jet reconstruction, indicating the influ- 
ence of the detector on the measurement. Results from 
a GEANT-based STAR detector simulation show good 
agreement with the data, demonstrating the reliability of 
the simulation framework used in the present analysis. 

Figure |2jb) shows the distribution of z for one of the 
bins in pion px in comparison to pythia with a GEANT- 
based detector simulation. To maximize the statistics in 
the simulation, the generator-level 7r°'s were used with- 
out requiring an explicit reconstruction. This led to a 
softening of the falling edge of the distribution at high 
z in simulations, since a full GEANT simulation was used 
for the containing jets, but did not affect the mean of the 
distribution. A small fraction of the events had z > 1, 
apparently corresponding to pions that carried more trans- 
verse momentum than their containing jet. This excess 
was caused by corrections applied during jet reconstruc- 
tion, which in some cases led to an underestimation of 
the jet energy, and was well reproduced in simulations. 
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FIG. 3: (color online) Longitudinal double-spin asymmetry 
for inclusive 7r° production at midrapidity in p + p collisions at 
y/s = 200 GeV, compared to NLO pQCD calculations based 
on the gluon distributions from the GRSV [33], GS-C [34] , 
and DSSV 12 global analyses. The systematic error (shaded 
band) does not include a 9.4% normalization uncertainty due 
to the beam polarization measurement. 



The asymmetry [Eq. ([TJ] was calculated as 

A = 1 ( N++ ~ RN+ ~) ™ 
LL P X P 2 (N++ + R.N+-) 1 K > 

where N ++ and A H are the tt° yields in equal and 
opposite beam helicity configurations, respectively, and 
R is the luminosity ratio for those two helicities. Typical 
values of R, measured with the BBCs to a statistical 
precision of 10 -3 -10 -4 per run, ranged from 0.85 to 1.2, 
depending on fill and bunch pattern. Figure [3] shows 
the measured longitudinal double-spin asymmetry for ir° 
production. The data points are plotted at the mean pion 
Pt in each bin. The lowest-pr point at 4.17 GeV/c was 
obtained from HT1 triggers only; other points are the 
HT1 and HT2 combined results. 

The systematic errors shown in the figure include 
point-to-point contributions from ir° yield extraction 
[(4— 14)xl0~ 3 ], invariant mass background subtraction 
[(6-11) xlO -3 ], and remaining beam background [(1- 
9)xl0 -3 ], as well as p^-correlated contributions from 
relative luminosity uncertainties (9 x 10~ 4 ) and from non- 
longitudinal spin components (3 x 10~ 4 ). All of the errors 
above are absolute errors on the measured asymmetry. 
An evaluation of the effects of non-longitudinal compo- 
nents of the beam polarization was not possible due to 
the limited statistics of 7r°'s in data taken with trans- 
versely polarized beams. Instead, the largest value from 
the jet measurement |llj over the relevant momentum 
range was taken as an estimate of this systematic error. 
An overall normalization uncertainty of 9.4% due to the 
uncertainty in the RHIC CNI polarimeter calibration is 
not shown. Studies of parity-violating single spin asym- 
metries and randomized spin patterns showed no evidence 
of bunch-to-bunch or fill-to- fill systematics. 
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In Fig. [3J the measured values for All are compared 
to NLO pQCD calculations [8] based on various sets 
of polarized gluon distribution functions. The DSSV 
curve [T2] is the result of the first global analysis that 
includes semi-inclusive and inclusive DIS data, as well 
as results obtained by the PHENIX [10] and STAR [TT] 
experiments. The GS-C curve [53] refers to a polarized 
gluon distribution function that has a large positive gluon 
polarization at low x, a node near x ~ 0.1, and a nega- 
tive gluon polarization at large x. The GRSV standard 
curve is based on the best fit to DIS data [55] . while the 
other GRSV curves show scenarios of extreme positive 
(Ag = +<?), extreme negative (Ag = —g), and vanishing 
(Ag = 0) gluon polarization at the starting scale [9] [33] . A 
maximal positive gluon polarization scenario, which has a 
total gluon spin contribution AG = J^Ag(x) dx — 1.26 at 

an initial scale of 0.4 GeV 2 [33J[35], is excluded by our mea- 
surement at 98% confidence level, including systematic 
uncertainties. This is in agreement with the conclusions 
from the inclusive jet measurements by STAR [TTJ |2"B] 
and from the inclusive ir° measurement by PHENIX [TP] , 
The data are consistent with all other gluon polarization 
scenarios, in particular with the DSSV case. 

In summary, we report a measurement of the invariant 
cross section and the longitudinal double-spin asymme- 
try All for inclusive tt° production at midrapidity with 
the STAR detector at RHIC. The cross section was 
determined for 1 < px < 17 GeV/c and found to be in 
agreement with a NLO pQCD calculation based on the 
CTEQ6M parton distribution functions and the DSS 
fragmentation functions. This set of fragmentation func- 
tions was constrained by data that included measurements 
of 7T° production at midrapidity by PHENIX JTU] and 
at forward rapidity by STAR |27j . The mean transverse 
momentum fraction of 7r°'s in electromagnetically trig- 



gered jets was found to be approximately 0.7 and to 
rise slightly with pion px, in agreement with a PYTHIA- 
based Monte Carlo simulation that included detector 
effects. This measurement has the potential to contribute 
to future fragmentation function studies. The asymme- 
try All was measured in the hard scattering regime at 
3.7 < pt < 11 GeV/c and found to be consistent with 
NLO pQCD calculations utilizing polarized quark and 
gluon distributions from inclusive and semi-inclusive DIS 
data and from polarized proton data. Our data exclude 
a maximal positive gluon polarization in the nuclcon, in 
agreement with results obtained from inclusive jet pro- 
duction in polarized proton collisions by STAR [TT1 12"8"] . 
while being a statistically independent measurement, sub- 
ject to a different set of systematic uncertainties. With 
increasing integrated luminosity, the neutral pion channel 
has the potential to provide additional constraints on the 
gluon polarization in the polarized proton. 
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